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The hydrodynamics, describing dynamical effects in superfluid neutron stars, essentially differs 
from the standard one-fluid hydrodynamics. In particular, we have four bulk viscosity coefficients in 
the theory instead of one. In this paper we calculate these coefficients, for the first time, assuming 
they are due to non-equilibrium beta-processes (such as modified or direct Urea process). The results 
of our analysis are used to estimate characteristic damping times of sound waves in superfluid neutron 
stars. It is demonstrated that all four bulk viscosity coefficients lead to comparable dissipation of 
sound waves and should be considered on the same footing. 
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I. INTRODUCTION 

The matter in pulsating neutron stars is not (even lo- 
cally) in chemical equilibrium. Particles of different kinds 
turn into one another so that the system evolves to equi- 
librium. If a deviation from the equilibrium is small then 
the processes of mutual transformations of particles can 
be described in terms of an effective bulk viscosity (see, 
e.g., Ref. [![). This viscosity influences the 'instability 
windows', that are the regions of physical parameters 
(e.g., the rotation period and temperature of a star) at 
which the neutron star becomes unstable against the ra- 
diation of gravitational waves 0,1 ill- The bulk vis- 
cosity, generated by non-equilibrium processes of particle 
transformations, was calculated in a series of papers for 
neutron-star matter of various composition (for example, 
for matter composed of neutrons, protons, and electrons 
with an admixture of muons; for hyperon or quark mat- 
ter). A short review and references to these papers can 
be found in Ref. @. 

It is generally agreed that the stellar matter becomes 
superfluid at a certain stage of neutron star thermal evo- 
lution 0, H, Q . A lot of attention has been paid to the 
question how superfluidity affects the bulk viscosity (see, 
e.g., Refs. 0, El, 03, OS QUI). In analogy with the or- 
dinary hydrodynamics of non-superfluid liquid, the only 
one 'standard' bulk viscosity coefficient has been calcu- 
lated and analyzed in all these papers. Meanwhile, it is 
well known [H. Il6l.[r^ that a superfluid liquid, composed 
of identical particles, is generally described by the four 
bulk viscosity coefficients. So, what can be expected from 
neutron stars, which contain a mixture of many super- 
fluid species? 

In this paper we show that non-equilibrium processes 
of particle transformations lead to the appearance of at 
least four bulk viscosity coefficients. Each of them is im- 
portant for analyzing dissipative processes in superfluid 
neutron stars. To be specific, we consider the simplest 
model of stellar matter composed of neutrons, protons, 
and electrons (npc-matter). In this case the bulk vis- 
cosity is associated with the non-equilibrium direct or 
modified Urea process. 

The paper is organized as follows. In Sec. [TT] we phe- 



nomenologically obtain the general form of the dissipa- 
tive corrections to the equations of relativistic hydrody- 
namics [l8l . [l9j , describing a superfluid liquid composed 
of identical particles. In Sec. IIIII this dissipative hydro- 
dynamics is generalized to describe superfluid mixtures 
and applied to npe-matter. In Sec. IIVI we calculate and 
analyse all four bulk viscosity coefficients provided by 
non-equilibrium beta-processes. For illustration of these 
results, in Sec.|V] we calculate the characteristic damping 
times of sound waves in superfluid npe-matter. Sec. I VII 
presents summary. 



II. THE DISSIPATIVE RELATIVISTIC 
HYDRODYNAMICS OF ONE-COMPONENT 
SUPERFLUID LIQUID 

In this section we obtain the general form of dissipa- 
tive terms entering the equations of relativistic superfluid 
hydrodynamics of electrically neutral liquid composed of 
identical particles. For that purpose, we need to choose 
a version of non-dissipative hydrodynamics. There is a 
number of equivalent formulations of non-dissipative rel- 
ativistic superfluid hydrodynamics [H HH, [2I, [H, Hi] . 
The most elegant and general am ong st them seems to be 
the formulation of Carter [H HH, l22| in which the hydro- 
dynamic equations follow from a convective variational 
principle. Most of the relativistic calculations (see Refs. 
[25l [2a , I27I ]) modelling pulsations of superfluid neutron 
stars have been made within this approach (see, however, 
0). 

In this paper, we will not use the Carter's hydrody- 
namics because it is an essentially phenomenological the- 
ory and it does not allow easy interpretation in terms of 
quantities calculated from microscopic theory. Since our 
main goal is the calculation of bulk viscosity coefficients, 
we will employ the hydrodynamics of Son [l8l |. It was 
initially proposed in the context of heavy-ion collisions 
and is derived directly from microscopic theory. There- 
fore, it is straightforward to relate various parameters of 
this hydrodynamics to microphysics. Using the notations 
of Ref. [n| it can be rewritten in a particularly simple 
form, which is a natural relativistic generalization of the 
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standard non-relativistic superfluid hydrodynamics pio- 
neered by Tisza [28| , Landau [2^, [3(| , and Khalatnikov 

a 

Although the Son's description is ideal for comparing 
with microphysics, it has some serious disadvantages. In 
contrast to the Carter's hydrodynamics, in which the ba- 
sic fluid variables are the particle number density current 
and the entropy density current, the Son's hydrodynam- 
ics is a hybrid in a sense that its fluid variables are the 
rescaled entropy density current and the rescaled momen- 
tum of a particle (or a Cooper pair) from the condensate 
(in the literature it is traditionally and somewhat confus- 
edly referred to as 'the superfluid velocity'). As a con- 
sequence, the Landau-type hydrodynamics of Son has a 
lower symmetry than that of Carter. However, it can be 
simply translated into the hydrodynamics of Carter as 
was demonstrated, for example, in the review paper by 
Andersson and Comer [33] (see their section 16.2). 

Below we will use the hydrodynamics of Son [18], em- 
ploying the notations of Ref. [191 ] . convenient for our 
problem. 

Unless otherwise stated, the speed of light is set equal 
to c = 1. 

The hydrodynamic equations take the standard form 



= 0, 



(1) 

(2) 



where 9 M = d/dx^ (space-time indices are denoted by 
Greek letters). Neglecting dissipation, the particle cur- 
rent density and the energy-momentum tensor T^ v 
can be presented in the form [191 ] 

f = nu^ + Yw^, (3) 
+Y {w^w v + fj, w^u v + nw v u»). (4) 



limit spatial components of the four-vectors u M and 
are equal to 



u = V 



w 



m(V s -V q ), 



(7) 



where V s = V<p/m and Vq are, respectively, the su- 
perfluid and normal velocities of the well known non- 
relativistic theory of superfluid liquids (see, e.g., Ref. 
0). 

The four-velocity must satisfy one additional con- 
straint. This constraint fixes a comoving frame, that is 
the frame where we measure (and define) all the thermo- 
dynamic quantities. Choosing the constraint in the form 



= 0, 



(8) 



from Eqs. and ([U one obtains that in this particu- 
lar case comoving is the frame where four-velocity equals 
•vP = (1, 0, 0, 0). It is straightforward to show that in this 
frame the basic thermodynamic quantities n, V<p, and e 
are defined by 



3 = n, 
3 = YV<i 
T 00 = e. 



(9) 
(10) 

(11) 



Using some equation of state and taking into account the 
second law of thermodynamics (T is the temperature, S 
is the entropy density) 

d£ = TdS + fidn+ — d(w' i w IJ ,), (12) 

as well as the definition of the pressure 

P = -e + [m + TS, (13) 



Here n, P, e, and \i (not to be confused with the space- 
time index fx\) are the number density, pressure, energy 
density, and relativistic chemical potential of particles, 
respectively; ip v = diag(— 1, +1, +1, +1) is the special- 
relativistic metric; Y is the relativistic analogue of su- 
perfluid density p s . In the non-relativistic limit, we have 
Y = p s /m 2 , where m is the mass of a free particle. Fur- 
ther, m m is the four-velocity of normal (non-superfluid) 
liquid component, normalized so that 

ttpu" = -1; (5) 

is the four-velocity, which describes the motion of 
superfluid liquid component. It can be expressed through 
some scalar function <p, 

up = &*<!> -nit". (6) 

It is easy to verify that <fi is related to the wave function 
phase $ of the Cooper-pair condensate by the equality 
= (H/2) V$ (see Ref. [HI). In the non-relativistic 



we can express all other thermodynamic quantities as 
functions of e, n, and V0. Eqs. (|T|)-(|13p fully describe the 
non-dissipative relativistic hydrodynamics of uncharged 
superfluid liquid. As a consequence of these formulae, one 
can easily derive the continuity equation for the entropy, 

fyOSu") = 0. (14) 

Now let us include dissipation in the hydrodynam- 
ics described above. As in the non-dissipative case, 
we assume that in the comoving frame [in which 
= (1,0,0,0)], the basic thermodynamic quantities 
e, n, and V</> are still defined by Eqs. (|9"))- (fTTj) . Being 
written in a relativistically invariant form, the conditions 
<j9j) and (fTQ|) imply that the particle current density j M is 
still given by Eq. ([3]), where some four- velocity w 11 sat- 
isfies the constraint © [as in the non-dissipative case]. 
In view of Eqs. 1(5)) and (fTUjl . the most general expression 
for the four-velocity is 

up = d^cj) - (ji + x)u". (15) 
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Here a scalar x is a small dissipative correction to be 
determined below. If we neglect dissipation, then k = 
and Eq. (fT"5|) coincides naturally with ([6]). 

Without any loss of generality, the expression for the 
energy-momentum tensor can be presented in the 
form: 

T"" = (P + e) u"u u + Pif v 

+Y{w fI w u + fiiiAt" +fjtw v u'')+T> iV . (16) 

Here t^ v is an unknown dissipative tensor. In view of 
Eqs. ([5]) and (fTTj) , t^ v satisfies the constraint 

u^u v t^ v = 0. (17) 

Let us determine the dissipative corrections r MI/ and h 
assuming they are linear in small gradients of hydrody- 
namic variables. For this aim we need to derive an en- 
tropy generation equation. It can be easily obtained from 
Eqs. (0-© if we make use of Eqs. ©, ©, JTJJ), CEO), 
(USD, and d]), 

d»S» = -^(W)-t^(^) 

+Yw» £ d^T + u- Yw» - 0„ V (18) 

Here the entropy current density is given by 

S" = SV - ^ - J Yw", (19) 

and satisfies the natural constraint u^S^ = —S. 

Let us analyze the last two terms in Eq. (fT5|) . In ad- 
dition to a quadratic dependence on small gradients of 
hydrodynamic variables, they also depend on the four- 
velocity As follows from Eq. (J7J), the spatial part of 
the four-vector is proportional to the difference be- 
tween the superfluid and normal velocities. In a large 
variety of problems this difference is small (as a conse- 
quence, the time component w° is also small because of 
the constraint [5]). In particular, it cannot exceed some 
(not very large) critical value AVer at which superfluidity 
breaks down (see Sec. lIV[) . For instance, if we study small 
perturbations of matter which is initially at rest (in ther- 
modynamic equilibrium with i// 1 = and = 0), then 
the last two terms in Eq. (fT8"|) are much smaller than the 
first two terms (a typical example of such situation is pro- 
vided by sound waves in superfluid npe-matter, see Sec. 
ly]) . Below, we will neglect the last two terms in Eq. (fT8|) 
when obtaining the dissipative corrections t^ u and x. 

Moreover, in the dissipative corrections we will also 
neglect small dissipative terms, explicitly depending on 
. For example, we will neglect the terms of the form 
w^d^T and u^w^d^u 1 in the expressions for x and t^ u , 
respectively. An inclusion of these small terms would 
result in 13 kinetic coefficients describing dissipation in 
superfluid liquid. In the non-relativistic case the same 
approximation is used, for instance, in the textbook by 
Landau and Lifshitz Pi (see their §140) and in the mono- 
graph by Khalatnikov [17| ■ 



Since the entropy does not decrease, the right-hand 
side of Eq. (fT5)) must be positive. This requirement puts 
certain restrictions on a general form of t^ v and k (linear 
in gradients). The standard consideration (see, e.g., Ref. 
[33| ) shows that 

t"" = — k (H^ u v + H u "< u") (<9 7 T + Tu s dsu-y) 

- ^1 H*" <9 7 (Yw r ) - 6 d jU ^, (20) 
x = -bd^iYw^-Ud^. (21) 

In these equations = rf- v + u^u v \ k and rj are, re- 
spectively, the thermal conductivity and shear viscosity 
coefficients; £1,. ..,£4 are the bulk viscosity coefficients. 
From the Onsager symmetry principle it follows that 

£1 = U- (22) 

For the positive definitcness of the quadratic form in the 
right-hand side of Eq. (fT8|) it is necessary to have the 
kinetic coefficients k, 77, £2, and £3 positive and the coef- 
ficient £1 satisfying the inequality 

£i 2 < 66. (23) 

In the non-relativistic limit the dissipative hydrody- 
namics proposed here coincides with the well known the- 
ory of Khalatnikov (see, e.g., Refs. [3, E3, ISS] ) - For illus- 
tration, let us indicate how the bulk viscosity coefficients 
Cichir ■ -j£Kh4 of Khalatnikov are related to those intro- 
duced in this paper. It is easy to demonstrate that 

£jChl = — , £Kh2 = f2, (24) 

m 

£Kh3 = 6<h4 = — ■ (25) 

Thus, in this section we have constructed the relativistic 
dissipative hydrodynamics of a superfluid liquid, com- 
posed of identical particles. It should be noted, that the 
dissipation was first included into the hydrodynamics [Lsj 
by Pujol and Davesne [35]. However, it is difficult to 
use their dissipative hydrodynamics in applications. The 
point is that the authors do not specify the comoving 
frame, where they define thermodynamic quantities. It 
is easy to verify that the frame which is defined as co- 
moving in our paper, cannot serve as comoving in Ref. 




III. VISCOSITY IN SUPERFLUID MIXTURES 

Let us apply the general formulae obtained in Sec. |TT] 
to neutron star matter. As already mentioned in Sec. 
HI we consider the simplest model of neutron star cores 
composed of neutrons (n), protons (p), and electrons (e). 
All results of this and following sections can be easily 
generalized to the case of matter with more complicated 
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composition (e.g., npe-matter with admixture of muons 
or hyperon matter). 

It is generally agreed that as a neutron star cools down 
the neutrons and protons become superfluid in its core. 
In such a system we have three velocity fields (instead 
of two, as in the previous section). They are superfluid 
velocity of neutrons, superfluid velocity of protons, and 
normal velocity u M of neutron and proton Bogoliubov 
quasiparticles and electrons. In this section we do not 
consider the dissipative effects related to the diffusion of 
particles. Neglecting the diffusion, nuclcon Bogoliubov 
excitations and electrons move with the same velocity u^. 
In analogy with Sec. [Til it is convenient to introduce four- 



velocities 



and 



j/ n \ aim u^/ p s instead of superfluid velocities of 
neutrons and protons, respectively. The Son's version of 
non-dissipative hydrodynamics was extended to the case 
of npe-mixture in Re f. |19f (for earlier formulations see, 

e.g., Refs. UmiEliallEi). The main goal of this 

section is to include the viscous dissipative terms into 
the hydrodynamics (l9j . Below the subscripts i and k 
refer to nucleons, i,k — n, p. Unless otherwise stated, 
the summation is assumed over repeated nucleon indices 
i and k. 

The full set of hydrodynamic equations describing su- 
perfluid mixtures consists of i) energy-momentum con- 
servation law @ with the energy-momentum tensor T M " 
given by 



-Y, 



ik 



w (i) w (k) 



■IM w^ k) u v 



'■ (26) 



ii) particle conservation laws written for neutrons, pro- 
tons, and electrons (/ = n, p, e), 

Mo =0 ' 3{ i) =niui' + Y ik w{ k) , f {e) =n e u»; (27) 

iii) constraints on the four- velocities , 

= 0, (28) 
and iv) the second law of thermodynamics, 



de = T dS + \ii drii + fi c dn c H — — d 



(29) 



To take into account potentiality of superfluid motion, 
four-velocities should be expressed through some 
scalar functions <pi and presented in the form 



qiA" - (jm + Hiju 1, 



(30) 



Note that one can avoid introduction of these new func- 
tions <pi in the hydrodynamics of superfluid mixtures 
if one formulates the potentiality condition ([3"0]) in the 
equivalent way 



(31) 



Below we will use the latter formulation because it is 
more suitable for our purpose. In this approach four- 
velocities wfy should be treated as independent hydro- 
dynamic variables. 

In Eqs. (f2"6"]) - (|3"Tj) m and n\ are, respectively, the rel- 
ativistic chemical potential and the number density of 
particle species I = n, p, c; A^ is the four-potential of 
the electromagnetic field; qi is the electric charge of nu- 
cleon species i. Furthermore, Yi k = Y k i is a 2 x 2 sym- 
metric matrix which naturally appears in the theory as 
a generalization of the superfluid density to the case of 
superfluid mixtures. In the non-relativistic limit this ma- 
trix is related to the entrainment matrix pi k (see Refs. 
0, ISIil) IH) b y tne equality Y lk = p lk /(m t m k ), where 
mi is the mass of nucleon species i. The pressure P is 
defined in the same way as for non-superfluid npe-matter 
(compare with Eq. [Tj 



P 



-e + Hifii + fi e n e + TS. 



(32) 



The dissipative hydrodynamics formulated above differs 
from the hydrodynamics [19j , describing superfluid mix- 
tures, only by the dissipative terms t^ v , x n , and >r p . 
The general form of these terms can be found from the 
entropy generation equation which is analogous to Eq. 
([15]) [see Sec. HI], 



T 



+Y ik w^ k) ^ df+T + u" Y ik w^ k) ^ dviif,, (33) 



where the entropy density current is 
u,, ,,,, H; 



S^ 1 = Su» 



rp r ^ rp Y ik W^ k y 



(34) 



Using the requirement that the entropy does not de- 
crease, one can easily obtain the dissipative terms t^ u , 
x n , and k p from Eq. (|33p . 

= -k (H^ u v + H v ~< it") (<9 7 T + Tu s dgu-y) 



Hm H u6 I g + g 



ZuH» v d 1 Y lk w~ [k) 



= —%5i dfj 



YikW^ 



(35) 
(36) 
(37) 



Here, as in Sec. [TT] we omit small dissipative terms, ex- 
plicitly depending on and, in addition, we neglect 
particle diffusion. An inclusion of these terms would 
result in 19 kinetic coefficients as it has been recently 
demonstrated by Andersson and Comer [4(| for the case 
of npe-matter (they obtained their result using a non- 
relativistic version of the Carter's hydrodynamics). 

In Eqs. ([35])-([37]) £u, £2, %3i, £«, and £ 5i are the bulk 
viscosity coefficients (i = n, p). Some of them are related 
by the Onsager symmetry principle, 



fit — £,4i, ?3p — ^5n ■ 



(38) 
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In addition, for positive definiteness of the quadratic form 
in the right-hand side of Eq. ([55)) one needs the following 
inequalities 



K > 0, V > 0, &n > 0, 
2£ln£lp£3p + ?2S3n?5p 



£lp£3n 



Equations (j33|) -(f3"9 |) are derived under the assump- 
tion that electrons and protons can move independently. 
However, this is not the case since they are charged. Any 
macroscopic motion of electrons is accompanied by that 
of protons to ensure quasineutrality condition (see, e.g., 

n c = n p . (40) 



One can obtain then from the continuity equations J27J) 
for protons and electrons (neglecting small 'diffusive' 
terms), 



<9 M Y pk w^ k) =0. 



(41) 



In principle (if we are not interested in the distribution 
of the electromagnetic field, which couples together elec- 
trons and protons), we can use this equation instead of 
the constraints ([28| and (|3~lj) for protons. 

Equations ([55j) - (j59")) should be modified to take into 
account the conditions (|40|) and (|41|) . As follows from Eq. 
(|4"Tj) . the only bulk viscosity coefficients which contribute 
to and K n (see Eqs. [551 and 155)) are £i n , £2, &3n, and 
^4 n . Since we neglect the last two terms in the right- 
hand side of the entropy generation equation (|33p . only 
these four coefficients are responsible for dissipation of 
mechanical energy of macroscopic motion in superfluid 
npe-matter. 



IV. NON-EQUILIBRIUM BETA-PROCESSES 
AND THE CALCULATION OF BULK 
VISCOSITY 

In Sec. IIIII we phenomenologically considered npe- 
matter and found the viscous terms in the relativistic 
hydrodynamic equations for superfluid mixtures. When 
doing this, we have ignored the fact that because of non- 
equilibrium beta-processes, the number of particles in the 
system is not conserved. In this section we demonstrate 
that the effect of non-equilibrium beta-processes is equiv- 
alent to the appearance of four effective bulk viscosity 
coefficients £i n , £2, £311, and ^4 n in the hydrodynamic 
equations. 

In the npe-matter we have two types of beta-processes 
responsible for beta-equilibration. They are the direct 
Urea process and the modified Urea process. The power- 
ful direct Urea process is open only for some equations of 
state with large symmetry energy (and at large enough 
densities, when p ¥a <p Pp + p Fc , P^, Ppp, and p Fc being 



the Fermi momenta of neutrons, protons, and electrons, 
respectively). The direct and inverse reactions of this 
process have the form (see, e.g., Refs. 0,0, H|) 



& P > 0, 6 > 0, 

!p6 - &fc P > 0. (39) 



p 



(42) 



Here v and D stand for neutrino and antineutrino, re- 
spectively. If the direct Urea process is forbidden by 
momentum conservation, then the main mechanism of 
beta-equilibration is the modified Urea process 0, 0, 0] 

N + n^ N + p + e + v, N + p + e— >N + il + v. (43) 

An additional nucleon N = n, p here is needed to take an 
excess of momentum away and open the process. 

The full thermodynamic equilibrium includes beta- 
equilibrium and we have [44j 



Sp = fi n - fi p - fi c = 0. 



(44) 



In this case the number of direct reactions in a matter ele- 
ment per unit time is equal to the number of inverse reac- 
tions. Thus, the total number of particles of any species 
remains constant. In particular, the electron generation 
rate Ar (that is the net number of electrons generated 
in beta-reactions in a unit volume per unit time) is zero, 
Ar = 0. 

If we perturb the system, the condition (j44|) will not 
necessarily hold (Sfj, ^ 0) and the direct and inverse reac- 
tions will not precisely compensate each other (Ar ^ 0). 
In this paper we assume that the deviation from the equi- 
librium is small. 5p <C kp,T. Then Ar can be presented 
in the form pH til QJ, E| HI 



Ar = A S/i. 



(45) 



Here A is a function of various thermodynamic quanti- 
ties defined in the equilibrium state (e.g., of the temper- 
ature and the particle number densities). For superfluid 
npe-matter this function was calculated by Haensel, Lev- 
enfish, and Yakovlev [HJGj]]. Note, that these authors 
neglected the dependence of Ar on the scalars w^w^a 
though (in principle) they can be non-zero in thermody- 
namic equilibrium. However, it seems that the results of 
Refs. [ljj, EH are accurate as long as (as an example, we 
take the case of i = k = n) 



— ~ 771 n (U sn - Vq) < k B T, (46) 

m n 

where V^ n is the superfluid velocity of neutrons. A nu- 
merical estimate of Eq. (|4"6"|) gives 



1 /2 

|AV n | = |F sn -F q |« 0.01c L^A 



(47) 



For clarity, we introduce the velocity of light c in this 
condition. On the other hand, as follows from the Landau 
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criterion (see, e.g., Ref. |47|), superfluidity of neutrons 
breaks down if |AV"„| > |AV" rr |, where 



|AV C 



An 



PFh 



0.0003 c 



Vl0 9 K; \n D 



??o 



1/3 



(48) 

Here A n is the energy gap in the neutron dispersion re- 
lation; T cn is the critical temperature of neutron super- 
fluidity onset; uq = 0.16 fm -3 is the number density of 
nucleons in saturated nuclear matter. It is worth not- 
ing that the criterion (|48[) is actually an upper limit on 
| AV cr | . In reality, a superfluid state can be destroyed 
at much lower |AV cr | due to the formation of vortices 
in superfluid matter (m [48|. Comparing Eqs. (|4"T|) and 
(|48f one can see that if the neutrons are superfluid then 
at not very low temperatures the condition (|47[) is always 
justified. 

To calculate the bulk viscosity coefficients we will use 
the non-dissipative hydrodynamics of superfluid mixtures 
[H (see also Sec. EH). The energy-momentum tensor 
r£iiv £ or guch a hydrodynamics is given by Eq. ([26]) , while 
the four-velocity satisfies the conditions (|28[) and 
Notice, that the dissipative components r A " y and 



Hi in Eqs. (J26J) and ((31]) should be taken zero, t» v = 
and >Ci = 0. Further, we assume that the quasineutrality 
condition (|40|) is fulfilled in both equilibrated and non- 
equilibrated matter. Using Eq. ([40]) . the second law of 
thermodynamics for mixtures (|29p can be rewritten in 
the form: 



de = T dS + /i n dnb — Sfj, dn c 



Y, 



ik 



Finally, let us assume that the four-velocities w 



(49) 
: in 



equilibrium. 

To take the non-equilibrium beta-processes into con- 
sideration it is necessary to add corresponding sources 
in the right-hand sides of the continuity equations for 
electrons, protons, and neutrons, 



3^ (n c u») = AT, 



= at, 
= -Ar. 



(50) 
(51) 

(52) 



When writing Eqs. (|50|) (|52p we bear in mind that every 
neutron decay is accompanied by the appearance of an 
electron and a proton (see the reactions H5] and |4"3"|) . 

Taking into account the quasineutrality condition (gD]) , 
one gets from Eqs. ([5U]) and ([5Tj) the equality ([IT]) . Thus, 
Eq. (|4"T]) remains the same as in the absence of beta- 
processes. It is more convenient to use the continuity 
equation for baryons instead of Eqs. (|5ip and ([S"!?]) . be- 
cause non-equilibrium beta-processes do not influence the 
total number of baryons per unit volume, = n n + n p . 
Summing together Eqs. ([5T|) and ([5"2]) and using Eq. ([4"T]) . 
one obtains 



n h u^ + Y nk w^ k) 



0. 



(53) 



Let us assume that npe-matter is slightly perturbed 
out of thermodynamic equilibrium so that deviations 
from the equilibrium are small and one can linearize 
the hydrodynamic equations. Below we will work in 
the comoving (at one particular moment) frame asso- 
ciated with some element of npe-matter. In such a 
frame = (1,0,0,0). We can assume further that 
perturbations in the comoving frame depend on time t 
as exp(iw c i), where lo c is the frequency of perturbation 
(measured in this frame). 

From the normalization condition © and Eq. ([2"8"]l it 
follows that 



(54) 
(55) 



w° {i) = 0, 

d t u° = 0, d t w° {i) = 0. 



Using these equalities, one gets from Eqs. ([50]) and ([53 

dtn e + div (n e u) = AT, 
d t n h + div [n h u + Y nk W( k) ] 



0. 



(56) 
(57) 



Here u and are the spatial components of four- 

vectors u M and w^js , respectively. The number densi- 
ties of electrons n c and baryons rtb can be presented as 
n e = n C Q + 5ti c , n^, = nt,o + c^b, where n C Q and «bo 
are the equilibrium number densities, while 5n c and Sn^ 
are small non-equilibrium terms depending on time as 
exp(iuj c t). Here and hereafter the thermodynamic quan- 
tities related to the equilibrium state will be denoted by 
the subscript '0'. Using these notations as well as formula 
([4"5*]) and linearizing Eqs. ([56"]) - ([FT]) , we get 



Sn 



1 



[A 6/i - n e0 div(u)] , 



(58) 



8n h = -- — {ribodiv(u) + div [Y~ nk W( fc) ] } . (59) 

IUJ C 

Notice, that the chemical potential disbalance <5/i in Eq. 
([58]) depends on Sn c and Sn^. Actually, <5/i can gener- 
ally be presented as a function of rib, n e , T, and the 



scalars w 



(the proton number density is equal to 



the electron one, see the quasineutrality condition 140]). 
One can neglect the temperature dependence o f 5a in a 
strongly degenerate npe-matter (see, e.g., Refs. fig. l50l| ) . 
Moreover, since the scalars w^^w^ are of the second 
order smallness, their contribution to 5fi is also negligi- 
ble (we recall that = in equilibrium). Expanding 
Sfiirib, n e ) in Taylor series in the vicinity of its equilib- 
rium value (which is zero, £/i(nbo, n c o) — 0, see Eq. I44|). 
one obtains in the first approximation 

d6/j,(n h0 ,n cQ ) cWju(n b o, «co) , 
dfi(n h , n c ) = dn h H dn a 



dn h0 



On, 



eO 



( 60 ) 

Analogous formulae can be written for perturbations of 
pressure SP = P(nb,n c ) — Pq, neutron chemical poten- 
tial Sfj, n = /i n (nb,n c ) — /i n o, and energy density Se = 



7 



e(n h ,n c ) - e , 

dP(n h0 ,n c0 ) dP{n h0 ,n cQ ) 
oP = - drib H 5 on c , (61) 



On 



].»() 



dnn(nbo,n e o) , 

d/Ltn = — 0?T.b 



<9n b0 
de(n h0 ,n c0 ) 
dn hQ 



dn e0 
dfi n (n hQ ,n cQ ) 
dn cQ 



Snv 



Sn c , (62) 
(63) 



In the last equation we have neglected the term of the 
form [de(nbo,n e o)/dn e o] 8n c . From the second law of 
thermodynamics (|49l we have cfe(ribo, n c o)/dn C Q = — <5/j. 
Therefore, this term is quadratically small and can be 
omitted. 

Using Eqs. ([Sg ll - M one finds 



, If. i • / \ . <9<5/z(n b0 ,n o0 ) . 

dn c = — < « n e o w c div(tt) H = nboAdiv(w) 

F ^ dn hQ 



d8^(n h Q,n e o 
dn h0 



■ A div [Y nk w {k) ] | 



(64) 



where F = u>'^ + iuj c \d5pL(n^,n c o)/dn C Q. For non- 
equilibrium Urca-processes and typical (for neutron 
stars) pulsation frequencies ui c ~ 10 3 — 10 4 s -1 , we have 
(see, e.g., Refs. dUEl) uj c > A \d8(j,(n h o,n e o)/dn e0 \. In 
this case Eq. can be simplified by keeping only terms 
linear in A. The result can be written as 



8n e — Sn c i + 5n e 2, 
where the first term equals 



dn c i = div(u) 



(65) 



(66) 



and describes compression and decompression of the pul- 
sating matter. This term remains the same even in the 
absence of non-equilibrium beta-processes. The second 
term 6n e 2 is due to non-equilibrium beta-processes 



Sn c2 = — r <j n h0 



dSfj,(n h0 ,x e0 ) 



dn h0 



div(u) 



dSfi(n h0 ,n e0 ) 
dn ha 



div [Y nk w 



(*)J 



(67) 



Notice, that in this formula the partial derivative 
d5fj,(ribo, x c o)/dn\>o is taken at constant value of x c q = 
n e o/nbo- When obtaining Eq. (|6T[) we used the identity 

9*(n b o,x e o) 9*(n bo ,n e o) , d^(n h0 ,n c0 ) 
"bo ~ = ™bo k 1- n c0 ■ 



dn h0 



dn h0 



dn c0 



where SS? is an arbitrary function of ribo and n e o- 

Our further strategy is as follows. We take the energy- 
momentum tensor T^ v for superfluid mixtures (with 
r A " / =0) from Eq. ([26]) and expand all the thermodynamic 
quantities (e.g., the pressure P and the energy density 
e), which determine this tensor, around their equilib- 
rium values. Restricting ourselves to linear perturbation 



terms, we obtain for the tensor T^ v (in the comoving 
frame) , 



T 00 = £ + Se, 
T°i = T j0 = faoYih w[ ky 

T^ m = (P Q + SP) 5 jm . 



(69) 



Here, the spatial indices j and m are equal to 1, 2, 3; the 
relativistic entrainment matrix Yn~ is taken in equilib- 
rium; SP and Se are given by Eqs. (|61|) and (|63p. respec- 
tively. 

Let us assume for a while that there are no non- 
equilibrium beta-processes in the matter. In this case 
Sn C 2 = (see Eq. IB"?} and the matter is reversibly pulsat- 
ing around the equilibrium. Then the mechanical energy 
is not dissipating, and the entropy is conserved. Thus 
it is obvious that dissipative are only those terms in the 
tensor T^ v which are directly related to Sn e 2- Writing 
out these terms in the form of a separate tensor t£ n., we 
have 



oo 

'bulk 

0j 
r bulk 



0. 



j0 
r bulk 



jm 

r bulk — 



dP(n h0 ,n c0 ) 
dn c0 



>n c 2 Ojm- 



(70) 



This tensor can be easily rewritten in an arbitrary frame 
if we take into account Eq. (|67[). 



M „ _ A dP{n h0 ,n c0 ) 



bulk 2 

uj c Oti c q 
f ddn{n h0 ,n c0 ) 
1 dn h0 



-«bo 



dSfj,{n h0 , x c0 ) 



dn h0 



(71) 



Comparing the tensor rt, ulk with the phenomenological 
dissipative tensor t^ v (see Eq. we find the expres- 
sions for the effective bulk viscosity coefficients £i n and 
£2, generated by non-equilibrium beta-processes 



\_ dP(n h0 ,n c0 ) dSfijnbo, n c0 ) 
,.,2 



A 



6 = o «b0 



dn c0 dn h0 
dP{n ha ,n c0 ) 9<5^(n bo ,a; o) 



On, 



c0 



On 



b0 



(72) 
(73) 



Let us do the same with the potentiality condition (T5T|) 
on the four- velocity of neutrons ^( n y As a result, we 
obtain in the comoving frame the dissipative component 
>r n , appearing because of non-equilibrium beta-processes, 



gMnQ^bQ 1 neo) 

<9n c0 



Sn, 



c2- 



(74) 



In a fully covariant form, this component is given by Eq. 
where the effective bulk viscosity coefficients ^ n and 
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6ta are 

^3n 



X_ d/j, n (n h o,n c0 ) dSfi(n b0 ,n c0 ) 
u>l dn c0 dn hQ 

d/J. n («bo , «oo) d8fi(n h0 , x c0 ) 



-~0 «bO 



dn h0 



(75) 
(76) 



In terms of the bulk viscosity formalism this means that 
one can vanish the quadratic form in the right-hand side 
of Eq. (|80[) by an appropriate choice of these velocities. 
This is possible only if the equality ([75]) is satisfied. 



and, in addition, the condition (|41[) is taken into account. 

Thus, we have calculated the four effective bulk viscos- 
ity coefficients £i n , £2, £311, and As will be shown in 
Sec. El each of them makes a comparable contribution to 
characteristic damping times of mechanical energy. No- 
tice, that only the coefficient £2 is usually analyzed in 
the literature devoted to non-equilibrium beta-processes 
in superfluid matter. The expression (1731) for £2 coincides 
with earlier results (see, e.g., Refs. [id. 

Not all of the coefficients ([72 ]) -(|73" ]) and ([75" ]) -([75 ]) are 
independent. The coefficients £i n and ^4 n are equal be- 
cause of the Onsager principle l|38p . This can be shown 
if one applies the following relation for npe-matter (see, 
e.g., Ref. Hi}). 



dP(nb ,x c o) 
dx cQ 



'bO 



dd n(n h0 ,x cQ ) 
dn h0 



(77) 



Furthermore, it is easy to verify, that instead of one of 
the inequalities ([39]) relating the coefficients £i n , £2, and 
£311, we have the strict equality 



(78) 



It is fulfilled only for those non-equilibrium processes, 
for which the expansion (|65[) is valid. Therefore, we have 
only two independent bulk viscosity coefficients. 

To prove Eq. (|78p it is instructive to consider the en- 
tropy generation equation. Neglecting all the dissipative 
processes except for the non-equilibrium beta-processes 
(e.g., neglecting thermal conductivity, diffusion, shear 
viscosity), one can obtain from the hydrodynamics dis- 
cussed in this section 



8fi AT = XSfi 2 . 



(79) 



Here we made use of Eq. ([45]) . Since we are interested 
only in terms linear in A, we can substitute 8n e \ for 8n e 
into Eq. (|60|) which determines Sfi. 

On the other hand, the entropy generation equation in 
terms of the effective bulk viscosities takes the form (see 

Eqs.E3 H35M37J!, and ED, 



+(£m + &n) 



(80) 



Let us compare the right-hand sides of Eqs. ([79"]) and 
([SO]). It follows from Eqs. O and (JTHJ) that for any given 
m m it is always possible to choose four-velocities wf^ in 
such a way, that <5/i = and the entropy generation rate 
vanishes (at some point and at some particular moment). 



DAMPING OF SOUND WAVES IN 
SUPERFLUID NPE-MATTER 



Let us illustrate the results of previous sections by 
calculating characteristic damping times of sound waves 
propagating in a homogeneous superfluid npe-matter. 
For simplicity, we consider only the damping due to the 
effective bulk viscosity. Neglecting dissipation, the sound 
modes of superfluid npe-matter have been thoroughly in- 
vestigated starting from the pioneering paper by Epstein 
[HTI ] in which he argued that there would be two typ es of 
sound modes in neutron stars (see, e.g., Refs. 
In particular, Gusakov and Andersson [li| were the first 
who considered in full relativity sound modes in npe- 
matter at finite temperatures. Here we closely follow 
their analysis. The pulsation equations (81) and (82) 
of Ref. [l!| can be used to describe sound waves tak- 
ing into account dissipation. Thus, there is no need to 
derive these equations from the hydrodynamics of super- 
fluid mixtures (Sees. IIIII and IIV[) once again. Instead, 
we will rewrite them using the notations adopted in our 
paper. The result is 



9t [(Po +S )U + flnoYnk W(k)] 
d t [MnOU + tfl( n )] = -VSfln. 



-V8P, 



(81) 
(82) 



The first equation is a consequence of the relativistic Eu- 
ler equation, which can be derived from Eq. |(3J) in a stan- 
dard way (see, e.g., Ref. [l[). The second equation follows 
from the condition ([3"Tj) written for neutrons. To fully 
define the system, Eqs. (fSTj) and ([52)1 should be supple- 
mented by the condition (|4T|) . Using Eqs. ([54]) and (|55[) . 
this condition can be presented in the form 



div [y p jfei0( fc )] = 0. 



(83) 



Now assuming that all the perturbations are plane waves 
proportional to exp(iwi — ikr), one obtains the following 
compatibility condition for Eqs. ([81]) - ([83| [s = ui/k is the 
velocity of sound in units of c] 



ys 



Cis 2 



C 2 + 8A = 0, 



(84) 
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where 
V = 
d = 
C 2 = 



/^nO (^nn^pp 
P 



-^np-^pn) 



MnO^bO 

Po 



(fh - 7i - 712/ ) + 72 - A> 
(/327i - A.72) , 



(85) 
(86) 
(87) 



nbo dP(n h0 ,x c0 ) _ n^o dfj, n (n ho ,x c o) 

7l 7-) ^ ) 72 ! 

rb C«bo Mno on h0 

n hQ dP(n h0 ,n e0 ) n h0 dfj, n (n h0 ,n e0 ) 

0i = — 75 , P2 = 75 ■ (89) 

P on h0 /ino on b0 

A small complex term SA appears in the compatibility 
condition (|84|) because of the bulk viscosity. It is given 
by 



SA 

A! 



s 2 A 2 ) 



MnO^bO 72 £ln - A*nO 02 £2 
--Po«bO 71 6n + -Po /3l ?4n, 



(90) 



(91) 

A 2 = UnO (y 6 + 6 + «b0 &n - ™b0 £ln - ™b0 Un) -(92) 

We remind that the bulk viscosity coefficients (and the 
quantities A\ and A 2 ) depend on the frequency to, Ai >2 ~ 
lo~ 2 . The biquadratic equation ([84)) has two non-trivial 
solutions for two possible sound velocities. Neglecting 
dissipation, these modes have been analyzed in details in 
Ref. [l9| . In particular, the sound velocities s± and s 2 ^ 
have been calculated there for the first and second modes. 
The dissipation leads to the appearance of small complex 
corrections Ssi, 2 to the velocities and consequently 
to decrements of sound waves. Since SA is small in com- 
parison with other terms in Eq. (|84|) . one can use the 
perturbation theory in deriving the characteristic damp- 
ing times ri ; 2. The parameters T\ and t 2 are e-folding 
times of the pulsation amplitude for the first and second 
sound modes, respectively, 



Tl 2 



k S 



81,2 



2 i s^ J 
' kSA 



(0)2 
2 



Ci 



(93) 



As follows from Eqs. (|9"01) - (f9"2"|) . they are independent of ui. 

At T — > T cn we have T nn ,Y np ,Yp n — > and y ks 
tbo/CMno^nn) — > oo (see Ref. 19] for a more detailed dis- 
cussion). In this limit the characteristic damping times 
are 



T2 



2Pq7i 
^2 



2 MnoPp Tl g 
w 2 (jiAi +fx a oD£ 2 ) 



(94) 
(95) 



Here we introduce the parameter D = 2 1i — 0il2- At 
T > T cn the neutrons are non-superfluid. In this case the 
second mode does not exist [formally, s 2 °^ = 0] , while the 
first mode is the usual sound wave. The characteristic 
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FIG. 1: Characteristic damping times Ti.,2 of sound waves 
versus temperature T for the first mode (two upper curves) 
and for the second mode (two lower curves) . The solid curves 
demonstrate the damping of sound taking into account all 
four bulk viscosity coefficients. The dashed curves are calcu- 
lated assuming that only £2 is non-zero. The neutron critical 
temperature is indicated by the vertical dot-dashed line. The 
baryon number density is — 3no = 0.48 fm~ 3 . 



damping time for an ordinary sound wave is given by 
Eq. (|94p. As expected, its damping is governed by the 
only one bulk viscosity coefficient £2- 

For illustration, in Fig. 1 we present the character- 
istic damping times n,2 of sound waves (in years) as a 
function of temperature T for two sound modes. The 
figure is plotted for npe-matter with the baryon number 
density nbo = Suq. The critical temperature of neutrons 
is taken to be T cn = 10 9 K. The protons are assumed 
to be non-superfluid. When calculating thermodynamic 
quantities and their derivatives we employed the equa- 
tion of state from Ref. [52j. It opens the direct Urea 
process at baryon number density of 5.84no. Therefore, 
the process is forbidden for ribo — 3no. In this case, 
the main mechanism of energy dissipation is the non- 
equilibrium modified Urea process. To calculate the func- 
tion A, which enters Eqs. l[?2" ]) -(j73" l) and l[75 ]) -(|7rj ]) for the 
bulk viscosity coefficients, we have used the results of Ref. 
11]. For calculating the relativistic entrainment matrix 
Yik we have employed the BJ v6 nucleon-nucleon poten- 
tial [Til [EH ■ Actually, the microphysics input we use 
here to plot the figure is taken from Ref. [lj| (see this 
reference for more details). 

The two upper curves correspond to the first sound 
mode, while the two lower curves - to the second mode. 
The solid curves are plotted taking into account all four 
bulk viscosity coefficients. The dashed curves are ob- 
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tained under the assumption that all coefficients but £2 
are equal to zero, £i n = ^ n = ^4 n = 0. 

As seen from the figure, the characteristic damping 
times of sound waves increase as the temperature de- 
creases. This is natural, because when the neutrons are 
superfluid, the Urea processes (and hence the bulk vis- 
cosity) are exponentially suppressed at T < T cn . Let 
us emphasize that at temperatures T < 5 x 10 8 K, the 
shear viscosity of electrons can exceed the bulk viscosity 
generated by the non-equilibrium modified Urea process. 
As a result, the damping of sound waves will be mainly 
due to the shear viscosity. 

The first mode turns into the ordinary sound at T > 
T cn . As follows from the figure, in the vicinity of neutron 
critical temperature the dissipation is primarily deter- 
mined by the bulk viscosity coefficient £2 (in accordance 
with Eq. I94p . Consequently, near the transition point the 
solid and the dashed curves for the first mode coincide. 
On the contrary, the difference between the solid and the 
dashed curves for the second mode remains significant at 
any T < T cn . The characteristic damping times for these 
two curves differ approximately by a factor of 3. 

It is worth noting that we would come to the similar 
conclusions if we considered sound waves in denser mat- 
ter, where the direct Urea process is open. In that case 
the characteristic damping times would be 6-7 orders of 
magnitude smaller, but the relative difference between 
the solid and the dashed curves will be approximately 
the same. 

Therefore, the main result of the present section is 
that the bulk viscosity coefficients £i n , ^3„, and ^4„ es- 
sentially influence the dissipative properties of superfluid 
npe-matter and cannot be ignored. All four bulk viscos- 
ity coefficients should be considered on the same footing. 

VI. SUMMARY 

We performed a self-consistent analysis of the influ- 
ence of non-equilibrium beta-processes on dissipation of 
mechanical energy in superfluid matter of neutron stars. 
We start with the Son's version of non-dissipative one- 
fluid relativistic hydrod yna mics to describe superfluid 
mixtures (see Refs. [18|, [l9|). We determined the gen- 
eral form of dissipative terms entering the equations of 
this hydrodynamics. For simplicity, the effects of parti- 
cle diffusion were ignored. The equations of dissipative 
hydrodynamics were applied to the matter composed of 
neutrons, protons, and electrons (npe-matter). In this 



case the hydrodynamic equations contain four bulk vis- 
cosity coefficients rather than one, as in non-supcrfluid 
matter. 

It was demonstrated, that non-equilibrium beta- 
processes generate all four bulk viscosity coefficients, and 
only two of them are independent. The other two coef- 
ficients can be expressed through the first two by Eqs. 
([55]) and ([75)1 . It is worth to emphasize that only the 
bulk viscosity coefficient £2 has been considered in the 
astrophysical literature so far. The expression ([75]) for £ 2 
coincides with similar expressions of previous works (see, 
e.g., Refs. [II [II 51, Si). 

To illustrate the results obtained in the present paper 
we considered a problem of damping of sound waves via 
the bulk viscosity due to non-equilibrium beta-processes 
in superfluid homogeneous npe-matter. It was shown 
that all four bulk viscosity coefficients make compara- 
ble contributions to the characteristic damping times of 
sound waves. 

Our results can be important for the analysis of vari- 
ous gravitational-driven instabilities in neutron stars, in 
particular, the r-mode instability (see, e.g., Ref. Q). The 
additional bulk viscosity coefficients lead to a more effec- 
tive damping of these instabilities. Moreover, the results 
can be applied to the problems of rotochemical and gravi- 
tochemical heating of millisecond pulsars with superfluid 
cores. In the absence of superfluidit y th ese problems were 
carefully analyzed in Refs. [H [Hill H|. The first at- 
tempt to discuss qualitatively the effects of superfluidity 
has been made in Ref. [57j . 

In conclusion let us note that the method of the bulk 
viscosity calculation, used here in the simple case of npe- 
matter, can be extended to matter with more compli- 
cated composition (npe-matter with admixture of muons, 
hyperon or quark matter). Such a generalization is be- 
yond the scope of the present study and will be consid- 
ered elsewhere. 
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